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Abstract 

The dispersive properties of non-degenerate optical parametric amplification in quasi-phase- 
matched (QPM) nonlinear quadratic crystals with an arbitrary grating profile are theoretically 
investigated in the no-pump-depletion limit. The spectral group delay curve of the amplifier is 
shown to be univocally determined by its spectral power gain curve through a Hilbert transform. 
Such a constraint has important implications on the propagation of spectrally-narrow optical pulses 
through the amplifier. In particular, it is shown that anomalous transit times, corresponding to 
superluminal or even negative group velocities, are possible near local minima of the spectral gain 
curve. A possible experimental observation of such effects using a QPM Lithium-Niobate crystal 
is suggested. 

PACS numbers: 42.65.Yj,42.65.ky,42.25.Bs 



1 



I. INTRODUCTION 



The theoretical and experimental study of second-order nonlinear processes has received 
in the past years a renewed interest after the introduction of the quasi-phase-matching 
(QPM) technique [I], which has lead to a major advance in applications such as frequency 
conversion, parametric oscillation and amplification, nonlinear optical frequency mixing and 
pulse shaping (see, for instance, 0, |j, § an d references therein). The QPM technique uses 
a periodic modulation of the nonlinear x coefficient (QPM grating) to compensate for 
refractive- index dispersion, permitting phase matching operation at any wavelength at room 
temperature, which would not be possible with birefringent phase matching. In addition, 
QPM enables the use of materials with strong nonlinearities which are not phase matchable 
by angle or temperature tuning, with a polarization state of interacting fields corresponding 
to the largest diagonal element of the x^ tensor. Though the physics of the QPM technique 
is known since a long time fTJ, only with the recent technological advances in the electric- 
field poling of ferroelectric materials, such as LiNbOs, LiTaOa and KTP, the experimental 
potentialities of QPM-based devices have become practicable and have reached nowadays a 
satisfactory degree of maturity p| . The possibility of engineering a QPM grating by breaking 
its periodicity introduces additional degrees of freedom for light control which open the way 
for a novel class of devices |], [7], [8], ||, [TO], [11], [L^, [13], [14], |15j. In particular, pulse compression 
and shaping by second harmonic generation (SHG) or by difference frequency generation 
(DFG) have been proposed and demonstrated using nonuniform QPM gratings || [7|, ||, [T2|, 
|I3"| . In such cases, it was shown that in the low-conversion regime the spectral transfer 
function governing the frequency conversion process is related to the longitudinal QPM 
grating profile by a simple Fourier transform (for a detailed analysis and for a discussion of 
the distinctive aspects of SHG and DFG cases, see [[12|, |13[ ) . The use of a QPM nonlinear 
crystal for parametric amplification of a signal field at frequency uj\ by a strong pump field at 



frequency has been demonstrated as well, both in waveguide and bulk geometries [16, 17 



with the achievements of gains as large as 40 dB using nanosecond pump pulses JT^]. As 
compared to most common parametric amplifiers based on birefringent phase matching, the 
use of a QPM-based parametric amplifier may offer the possibility of molding the spectral 
gain response of the amplifier by a suitable design of the QPM grating (see, e.g., [|l5|j). 
Though the process of parametric amplification is accompanied by the generation of the 



2 



idler field at frequency ui 2 = oos — u\ and is thus equivalent to the DFG scheme recently 



studied in literature |T3|, P^| , the simple Fourier analysis developed for DFG |I3]| , which 
neglects amplification of the signal field, is inadequate to study the parametric amplification 
process in QPM gratings, and a more complex analysis is needed, based e.g. on matrix 
transfer techniques [jn| . A general theoretical analysis of the properties of a QPM parametric 
amplifier, and especially of the interplay between its gain and dispersive properties, has not 
been pursued yet. Such a problem is closely related to the "inverse problem" of QPM grating 
synthesis, i.e. of the determination of a QPM grating profile that realizes a desired spectral 
gain response of the amplifier. 

In this paper we provide a general analysis of the process of parametric amplification in 
an aperiodic QPM nonlinear crystal. We will show in a very general way that the spectral 
power gain curve of the amplifier univocally determines its dispersive properties, whereas the 
QPM grating profile that realizes such a spectral gain curve is not uniquely determined. The 
interplay between the spectral power gain and dispersive curves of the amplifier may lead 
to abnormal dispersive properties of the amplifier, such as the occurrence of superluminal 
or even negative group velocities similar to those observed in atomic amplifiers with a gain 



doublet |18|, [T^, |20|. The paper is organized as follows. In Sec. II the basic equations 
describing parametric pulse amplification in a QPM nonlinear crystal are reviewed. In the 
undepleted pump approximation, general properties of spectral gain and dispersion curves 
are derived in Sec. Ill, and the inverse problem of QPM synthesis is addressed. Section 
IV deals with the propagation of spectrally narrow pulses and analyzes the occurrence of 
anomalous regions corresponding to superluminal or even negative group velocities. Finally, 
in Sec.V the main conclusions are outlined. 



II. PARAMETRIC AMPLIFICATION IN QPM QUADRATIC MEDIA: BASIC 
EQUATIONS 

The starting point of the analysis is provided by the propagation scalar wave equation 
for a linearly-polarized electric field £(z,t) in a nonlinear x medium with a quasiperiodic 
QPM grating profile. Under the usual plane-wave approximation and taking into account 
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material dispersion, from Maxwell's equations one can write (see, for instance, pl| ) 
d 2 £ r , ,„ N ~, . . _ , d 2 V NL 



dz 2 



/oo _ Q2'pl\L 

duk 2 (uj)£(z, u)) exp(-iujt) = /i ^ , (1) 



where 8(z,u) = (2n) 1 dco£(z, t) exp(iut) is the Fourier transform of S(z,t), k(uj) = 
(oj/cq) yl + x(oo) = {(jj/ca)n(oj) is the dispersion relation defined by the complex linear 



susceptibility x(uj) [or by the complex refractive index n(u>) = \J\ + x{ u )]> cq is the speed of 
light in vacuum, /xq is the vacuum magnetic permeability, and *p NL is the nonlinear driving 
polarization term. For a quadratic medium and neglecting dispersion effects of second- 
order polarization, one can take V NL (z,t) = eox ( z )E 2 (z, t), where x^ is the spatially- 
modulated nonlinear susceptibility that accounts for the QPM grating. To study parametric 
amplification or DFG, we assume that £{z,t) is described by the superposition of three 
wavetrains with carrier frequencies u)\, 0J2 and u 3 = Ui + U2, corresponding to the signal, 
idler and pump waves, respectively. We thus set 

£(z,t) = - [£i(z, t) exp(-iuit) +£2(2, t)exp(—iu 2 t) + £ 3 {z,t) exp(-iu 3 t) + c.c] , (2) 

where the envelopes £1,2,3(2, t) of wavetrains are assumed to vary slowly with respect to 
time t as compared to the exponential terms. In case of parametric amplification (or DFG), 
a strong pump wave and a weak signal wave are incident upon the nonlinear medium at 
the entrance plane z — 0, whereas the idler wave is generated by parametric interaction in 
the nonlinear medium and appears at the output plane z = L (see Fig.l). Substitution of 
Eq.(2) into Eq.(l) and setting equal the terms oscillating at the same frequency, yields the 
following set of coupled-wave equations: 

d 2 £\ , , 2 1 1 ' f\ \ c (2\ f u i X 



+ k\u x + id t )S 1 = -x (2) (^) £& (3a) 
ife * + k 2 (co 2 + id t )£ 2 = - X {2) ( ^) 2 £& (3b) 



dz 



OZ z \ Co 

d 2 £ 3 s-^ 2 



H + k\u 3 + id t )£ 3 = - X (2) (£) (3c) 

In deriving Eqs.(3), we have neglected the nonresonant terms in the nonlinear polarization 
driving term J> NL J and we used the following property (see, for instance, pif ): 



dwk 2 (uj)£(ijj) exp(— iut) = [k 2 (uo + id t )A(t)\ exp(— iu t) (4) 



4 



which is valid for any signal of the form £ it) = A(t) exp(-iuot), where the operator k 2 (u + 
idt) is defined through its power series expansion: 



k 2 (u + id t ) = k 2 (u ) + 2ik(u ) 



dk\ d 
du).. dt 



( dk\ 2 , . , / d 2 k 



,2 



U>0 



d 2 

gp + ... (5) 



Notice that, in case of a wavetrain at carrier frequency u and slowly- varying envelope A(t), 
the power series on the right hand side in Eq. (5) may be truncated at some order. At leading 
order in the expansion one retrieves the wavenumber at the carrier frequency Uq, the second 
one determines the group- velocity v g = \/{dk/doj) of the wavetrain, whereas higher order 
terms account for group velocity dispersion. If in Eqs.(3) the nonlinear terms are assumed 
to be weak and we truncate expansion (5) keeping the first term only, one sees that at 
leading order the amplitudes E\ are oscillating in space like exp(±ikiz), where k\ = k{uo{) 
(I = 1,2,3). If we consider only forward propagating waves, it is convenient to remove the 
fast oscillation over the wavelength spatial scale by setting Si(z,t) = Ai(z,t) exp(zA^z). For 
the QPM grating, we assume that is a quasi-periodic function of z with period A, i.e. : 

oo 

X (2) (z)= xi 2) Wexp(-2mvrz/A), (6) 

n=— oo 

where the Fourier coefficients Xn (z) are slowly varying functions of z over one period A. 
In practice, the slow dependence of coefficients on z can be achieved by a +/- reversal of 
domains in the ferroelectric crystal with a local period and local duty cycle that are slowly 
varying along the z axis (see, for instance, [12]). We assume that the nominal QPM period 



satisfies the phase matching condition: 



A = M§, (7) 



where Ak = k^ — k2 — k\ is the wavevector mismatch of interacting waves and M is a positive 
integer (QPM of order M). The evolution equations for the envelopes Ai(z,t) (I = 1,2,3) 
can be derived, in the limit of weak nonlinearity and quasi-monochromatic wavetrains, by a 
multiple-scales asymptotic expansion (see, for instance, [01); for the sake of completeness, 
a brief account of the derivation of the envelope equations is given in Appendix A. These 
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reads: 



f)Ai Ik 2 
2ik * ~jr = l k l ~ + id t)} A i - -^ d effA* 2 A 3 (8a) 

r?>4 2/r 2 
2^ 2 —^ = [A; 2 - A; 2 (cu 2 + idt)] A 2 - -^d eff A\A 3 (8b) 



dz L ' J n\ 

f)An 9k 2 
2*^ = fe 2 - A; 2 (cu 3 + A 3 - -l-dl H A x A 2 (8c) 



where: 



deff(z) = ^X (2) (^) exp(iAfe) = ^xS (9) 

and the overline denotes a spatial average over a few modulation periods of the QPM grat- 
ing. Equations (8) are the basic equations describing QPM parametric processes in the 
time domain. The linear operators on the right hand side of Eqs.(8), describing the linear 
dispersive and absorptive properties of the medium, may be expanded in power series [see 
Eq.(5)], and the number of terms that need to be kept depends on the the spectral exten- 
sion of the wavetrains. In the following, we will consider spectral regions of transparency 
for the medium, so that we will neglect the imaginary part of the dispersion relation k(u). 
In addition, for spectrally-narrow envelopes, only the leading-order term in the expansion 
may be kept, and Eqs.(8) reduce to the following ones: 

• 2d eff A* 2 A 3 (10a) 



dz 
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Vgl 


dt 


dA 2 
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dA 2 
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dt 


dA 3 
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dA 3 


dz 


Vg 3 


dt 



n 



= %%d eS! A\A 3 (10b) 
n 2 



k 



i- 



n 



d* eff A x A 2 (10c) 



where v g i = l/(dk/du) UJl (I = 1,2,3) are the group velocities of the three wavetrains. For 
spectrally broad wavetrains, higher-order terms in the expansion (5) should be considered, 
which are responsible for group velocity dispersion and higher-order dispersion effects. 



III. SPECTRAL GAIN CURVE AND GROUP DELAY ANALYSIS 

In this section we provide a general analytical framework to determine the complex spec- 
tral gain curve of a QPM parametric amplifier in the limit of undepleted pump and derive 
a general relation among dispersive and gain properties of the amplifier. We assume that 
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the amplifier is pumped by a continuous-wave (CW) or quasi CW strong pump wave at 
frequency us, and a weak signal field at carrier frequency u>i, described by the envelope 
Ai(0,t) = J™ f(Q) exp(— iQt), is incident upon the nonlinear crystal at the plane z — 0. 
In the limit of no pump depletion, which is valid for a low or moderate pump conversion, 
we may disregard Eq.(8c) and assume in Eqs.(8a) and (8b) a constant value of A 3 , given by 
A 3 = [2I S / (eoCons)] 1 ^ 2 , where I3 is the intensity of the incident CW pump and n 3 = n(cu 3 ). 
In this case, Eqs.(8a) and (8b) become linear equations, and the signal field envelope Ai(L, t) 
at the exit of the medium can be written using standard Fourier analysis in the form: 

/•oo 
dSlg(n)f(n) exp(-ifii) (11) 
-00 

where f(Q) is the Fourier transform of the incident signal wavetrain and g(Q) is the spectral 
gain curve of the amplifier. In order to determine the spectral gain curve g(fl), let us search 
for a solution of Eqs. (8a) and (8b) in the form: 

Ai(z,t) = u(z)exp[-ittt + i(3{VL)z\ (12a) 

A 2 (z,t) = v*(z) — J^exp[iQt-iP{Q)z], (12b) 
n 2 V m 

where Q is the frequency offset from the reference carrier frequency cu 1 of the signal wave, 
n\ = n(ui), n 2 = n(u 2 ), and (5 = f3(£l) is defined by the relation: 

m = w ± ^ 2{uJl + n) ~ k2 ^~w 2 ^ 2{uJ2 ~ n) ~ ^ ■ (13) 

The evolution equations of the complex envelopes u(z) and v(z) are obtained after inserting 
the Ansatz (12) into Eqs. (8). One then obtains for u(z) and v(z) the following coupled-mode 
equations: 

dn 

— = i5u + iq(z)v (14a) 
dz 

— = —i5v — iq*(z)u (14b) 
dz 

where 5 = S(Q) and q(z) are give by: 

m = w x ^ k2{uJl +Q) ~ k ^ + w 2 ^ 2{uJ2 ^ ' (15) 



q(z) = 2nd ef Az)J 2h xx . (16) 

" eoCbnira 2 n 3 AiA 2 V ; 
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In Eq.(16), Ai and A2 are the wavelengths in vacuum of signal and idler fields. Equations (14) 
have the form of the Zakharov-Shabat system encountered in problems of inverse scattering 
23| and grating analysis |2|, ^5], |26|, ^7|, the QPM grating profile d e ff(z) playing the role 



of the scattering potential [see Eq.(16)]. The general solution to Eqs.(14) can be written as 
(u(L, 8), v(L, 8)) T = M. (u(0, 8), v (0, 8)) T , where the elements of the 2x2 transfer matrix 
M = M{8) satisfy the conditions M22 = M* u , M21 = M* 12 , and detM = M n M 2 2 - 
M.yiM-21 = 1- From Eq.(12a), one realizes that the spectral gain curve g(fl) of the amplifier 
is given by [^8 : 



g (n) = M n (n) exp[iP(Q)L] (17) 

where the dependence of Ain on Q is through the function 8 = 8(Q) [see Eq.(15)]. The 
explicit calculation of M.(8) for a generic scattering potential q(z), i.e. QPM grating profile 
d e ff(z), can be done numerically by a standard cascading technique in which the grating is 
decomposed as the cascade of successive uniform sections and M. is derived as the product 
of the elementary matrices of each uniform grating section (see, for instance, [0). The 
application of the cascading technique for the study of DFG and parametric amplification in 
certain nonuniform QPM gratings has been recently addressed in [FJ, [Tj|]. The determination 



of the spectral gain curve of a QPM parametric amplifier for a given QPM profile is thus 
easy: first, one has to compute the matrix element M.\\ associated with the scattering 
problem [Eqs.(14)] by, e.g., the transfer matrix technique, for different values of frequency 
Q; the complex spectral gain curve is then retrieved by Eq.(17) with the use of Eq.(13). 
Here we are mainly concerned on investigating the interplay between amplitude and phase 
relationship of the complex spectral gain curve, which can be derived in a very general way 
using the results of the inverse scattering theory. To this aim, it is worth observing that, from 
inverse scattering theory, the scattering coefficient a(8) = M.n(5) exp(— i5L) is an analytic 
function in the lower part of the complex 5 plane, i.e. for Im(5) < 0, \a{8)\ > 1 on the real 
axis, and a(S) —>■ 1 as 5 —>■ 00 |p3| . In addition a(S) has no zeros for lm(5) < 0, so that 
l/a(S) belongs to the class of causal and minimal phase shift functions (see, for instance, 
|29f). This means that the knowledge of the modulus of a(5) univocally determines its phase. 
In fact, since the function F(S) = (da/ '88) /a = 8(\na)/85 is analytic in the lower part of 
the complex 5 plane and goes to zero for 8 — >• 00, the real and imaginary parts of F(8) can 
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be related by a Hilbert transform, that is: 

m = i r (18) 

After setting a(5) = \a(S)\ exp[i(f> a (6)] and making equal the imaginary parts of both sides 
in Eq.(18), one obtains: 

In order to further proceeds in the analysis, we must specify the dependence of (3 and 5 on 
f2. From Eqs.(13) and (15), the expansion of k — k(u) in in power series of f2 yields: 

^4K^ + i) + ^ 2 (^-^ + l-l) + --- (2o) 



5(Q) = ^ ( — - — ) + -tf ( k'l + k" 2 + % + f) + ... (21) 



2 \v gl v g2 J 4 \ - ki k 2 

where k' 12 and k'[ 2 are the first and second derivatives of k(uo) evaluated at the frequencies 
u)\ t 2- F° r the sake of simplicity, let us assume that v g \ ^ v g2 , which is usually the case for 
a non-degenerate interaction, and that the bandwidth of the amplifier is narrow enough to 
neglect group- velocity dispersion and higher-order dispersive effects in Eqs.(20) and (21). 
At leading order we may hence assume: 

« n > = Hi - £) • m ^ n (i + i) (22) 

If we introduce the spectral power gain curve G(Q) and group delay t 9 (Q) of the amplifier, 
defined as G(Q) = |g(f2)| 2 and r g (Q) = d(j) g /dQ, where <f> g (£!) is the phase of g(fl), from 
Eqs.(17), (19) and (22), and recalling that a(S) = M.u(8) exp(— i£L), one finally obtains the 
following relationship between the power spectral gain G and group delay r g of the amplifier: 



t,=—+\(- — -)- r ^-^W 1 ^ <23) 

v g i 2 \v gl v g2 J 7t J_ 00 85' 5' -5 

that is: 

^lnv/Gp 7 ) 1 

nz — i ui 

Vgl 7T J_ 00 

where the upper [lower] sign occurs if v g \ < v g 2 [v g \ > v g 2\- Equation (24) shows that the 
power spectral gain curve G{VL) of the amplifier uniquely determines its phase response, 
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T 9 



L x 1 r ^vGW) 1 m 



i.e. the group delay r g . The constraint imposed by Eq.(24) has important impacts on the 
response of the amplifier, reducing the flexibility in designing a given spectral gain curve as 
compared to the related DFG problem fl3|| . For instance, the possibility of "squaring" the 
spectral power gain curve of the amplifier comes at the price of an increased phase distortion 
. Such limitations come out just because the spectral gain curve satisfies the minimum- 



phase condition, which is not the case of the DFG transfer function [28 ] . Another interesting 
consequence of Eq. (24) is the existence of spectral regions where the group delay becomes 
superluminal or even negative, which will be studied in detail in the next section. Here we 
end our general analysis by briefly addressing the issue of the synthesis of the spectral gain 
curve g(fl) of a QPM amplifier (the so-called "inverse problem"), that is the determination 
of the QPM profile d e ff(z) that realizes a target spectral response g(Q). Such a problem has 
been extensively studied on several occasions in the context of optical filters and gratings 



theory (see, for instance, and references therein), and therefore it will be briefly quoted 
here. First of all, let us recall from inverse scattering theory that in Eq.(14) the scattering 
potential q(z) is uniquely determined from the knowledge of r(5) = M12/M.11, which must 
satisfy minimal requirements of causality, and that numerical techniques for determination 
of q(z) from the knowledge of r(5) are well developed |5I[ |32|. Let us now assign a target 
power spectral gain curve G(Q), and ask if there is a QPM grating profile d e ff(z) that 
realizes such a gain response. First of all, G(Q) should satisfy minimal conditions: G(Q) > 1 
for all real frequencies Q, be analytic in the lower [upper] part of the complex Q plane if 
v g l < v g 2 [v g l > V g ?\, G (ft) — > 1 for Q — > oo , and the integral J °° dfl\nG(Q)/ (1 + Q 2 ) should 
not diverge. In this case, the scattering coefficient a(5), and hence the matrix coefficient 
M 11, can be univocally determined owing to the minimal phase shift condition through 
Eq.(19). Once Ain(S) has been calculated from G(Q), the modulus of r(5) is determined 
by |r(5)| 2 = (1 — l/|7Wn| 2 ), however some degrees of freedom are left in the phase of r(5), 
yet ensuring the requirements of causality (see Sec. 4 of Ref . p9| ) . By changing the phase of 
r(8), one basically changes the spectral phase response of the generated idler wave, which 
is related to the phase of M.n (see J28J), without affecting the spectral gain amplitude and 
phase at the signal wavelength. For a physically realizable power spectral gain curve G(Q), 
there is not, hence, a unique QPM profile that realizes such a gain response, the different 
QPM profiles corresponding to different spectral phase response of the generated idler wave. 
The non- uniqueness of the inverse problem is analogous to that encountered, e.g., in the 



10 



design of fiber Bragg gratings for use in transmission |33|, |34 



IV. ANOMALOUS PROPAGATION OF SPECTRALLY NARROW OPTICAL 
PULSES: SUPERLUMINAL AND NEGATIVE TRANSIT TIMES 

A physically important consequence of the dispersive properties of a QPM parametric 
amplifier, stated by Eq.(24), is the existence of abnormal group delays, that is the existence 
of spectral regions where a spectrally-narrow signal pulse can travel through the amplifier 
with a superluminal {r g < L/cq) or even negative [r g < 0) transit time [j3"5[| . Let us consider 
the propagation of a spectrally-narrow signal pulse through the QPM amplifier. We assume 
that the spectrum f(fl) of the incident pulse is centered at a frequency Qq, so that we can 
write Ai(0,t) = exp(—iQ t)ho(t) and Ai(L,t) = exp(— iQ t)hL(t) at the input and output 
planes of the amplifier. If we now assume that the spectral extent of the envelope ho(t) is 
narrow enough such that the gain curve g(Q) of the amplifier varies slowly over the pulse 
bandwidth, an asymptotic expansion for hiif) may be obtained, as detailed in Appendix B. 
At leading order one obtains: 

h L (t) ~ g ho(t - At) (25) 

where go — fl'(^o) is the gain of the amplifier at the pulse central frequency Qo an d Ar is 
the complex group delay time, which is given by (see Appendix B): 



/ d\ng\ ,„ s . / d\n\/G 



ATS -'W , = T '^-^-gn-J 0j - (26) 

Notice that, since the group delay Ar is in general complex- valued, even at leading order 
in the expansion [see Eq. (B7) in Appendix B] the amplifier produces a pulse distortion p6 



However, if Q is a stationary point of the spectral gain curve, than the complex group 
delay Ar becomes real- valued and equal to the usual group delay r g = d(p g /dVL. In this 
case, at leading order the amplifier is distortionless and the pulse transit time is related to 
the spectral gain curve G(Q) by means of the Hilbert transform given in Eq.(24). Let us 
now suppose that v g \ > v 92 and that is a maximum of the spectral power gain curve; 
then for Q ~ Qq, the contribution to the integral on the right hand side in Eq.(24) near the 
singularity Q' = Q is negative, and the integral is thus expected to be negative. This implies 
that r g > L/v gi , i.e. the effective group velocity of the pulse is lower than v g \. Conversely, 
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if Qq is a minimum of the spectral power gain curve, the integral in Eq.(23) is expected to 
be positive, and the transit time r g lower than L/v gi . This means that superluminal or even 
negative transit times may be expected whenever the pulse spectrum is centered around a 
minimum of the spectral gain curve |37|| . The existence of anomalous transit times for a 
pulse tuned midway two amplification peaks is indeed a rather general result, and previous 
demonstrations of negative transit times have been reported for optical pulses as well as for 
electronic signals using a gain-doublet amplifier [|18|, [19|, |2(J, |38[ . 

To get a quantitative analysis of such effects, let us refer to a periodically-poled Lithium 
Niobate (PPLN) crystal pumped at the wavelength A3 = 532 nm with a signal field at 
Ai = 1.55 /im. The idler wave corresponds to the wavelength A2 = 810 nm; we further 
assume extraordinary wave propagation, so that d = d 33 ~ 27 pm/V. The temperature- 
dependent dispersion relation k = k(u) = n(u)u/co for extraordinary waves in Lithium 
Niobate is determined using Sellmeir equations from Ref . ||39|| . At 25°C, one can estimate 
v g i ~ 0.4815co, v 9 2 — 0.44220co, and a first-order QPM period A = 7.39 /im, which is 
accessible with current poling technology. The main task in the design of the QPM grating 
is to realize a gain curve G(Q) with a local minimum, e.g. at Qq = 0. There are several 
possibilities to achieve a spectral dip in the gain curve of a QPM amplifier, such as the 
inclusion of a defect in an otherwise periodic QPM grating or the cascading of two periodic 
QPM gratings. We consider here the latter case, however a similar analysis could be done 
for the former configuration. The QPM grating, shown in Fig. 2, consists of a sequence 
of two +/- square- wave uniform gratings, each of length a and period A = 2tt / 'Ak (first- 
order QPM), separated by a distance I. For such a structure one can easily calculate d e ff, 
and hence q(z), obtaining q(z) = q for < z < a, q(z) = for a < z < a + I, and 
q(z) = goexp(i$) for a + l < z < 2a + l, where go = (2/ K)d3 3 {8n 2 1 3 / ( y e conin2n 3 \i\2)] 1 / 2 and 
$ is a phase shift that depends on the relative phases of the two square waves in the two 
grating sections. For instance, the reversal of sign of d e /f, corresponding to $ = 7r, occurs if 
the two square waveforms in the two grating sections are shifted each other by half a period 
[see Fig. 2(b)]. Notice that, in the limiting case I — > 0, the structure reduces to a periodic 
QPM grating of length 2a with a defect at the center of the structure. In order to calculate 
the spectral gain curve g(£l) of the amplifier, we need to evaluate the matrix element Aiu 
according to Eq.(17). The transfer matrix of the structure shown in Fig. 2 can be readily 
calculated as the product of the three transfer matrices corresponding to the propagation 
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in the first uniform grating section, the middle grating-free section, and the second uniform 
grating section in the reverse order, i.e.: 

f exp(i<W) 



cosh (8a) + if sinh (8a) i exp(i<l>)^ sinh (0a) 
M = \ I x 

— i exp(— i$) 2& sinh (0a) cosh (8a) — if sinh (0a) 

cosh (8a) + if sinh (da) if sinh (0a) 

:<5 . 



^ exp(— iSl) 



x 



(27) 



— i^ 1 sinh (0a) cosh (da) — if sinh (8a) 

where ee ^M 2 - £ 2 - From Eqs.(17) and (27) one readily obtains: 

5 2 5 la I 2 

g(Q) = exp(iSl + if3L) cosh 2 (0a) - — sinh 2 (0a) + i- sinh(20a) + exp(i$)^- sinh 2 (0a) exp 

6 6 6 

(28) 

where (3 = P(Q) and 5 = 5(£l) are given, in the general case, by Eqs.(13) and (15). Around 
the resonance Q = 0, the approximate equations (22) may be used, which require the 
knowledge of the group velocities for signal and idler fields solely. A typical behavior of 
power spectral gain G(Q) and corresponding group delay r g (Q) for a few different values of 
$ are shown in Figs. 3,4 and 5. The parameter values of the QPM grating structure are a = 3 
mm and / = 2 mm. Notice that, in correspondence of local minima in the power spectral 
gain curve, the group delay curve shows local minima, corresponding to superluminal and 
even negative group delays (see, for instance, Fig. 3). An inspection of Eq.(28) reveals that, 
for $ = 7r, the curve G(£l) = \g(Q)\ 2 shows a dip at = 0, with G(0) = 1 and with a 
corresponding minimum in the group delay curve which can be calculated analytically in a 
closed form and reads: 



rW = 0) = L(±-±-)-'(±-±)x 



cosh 2 (|g |a) + sinh 2 (|g |a) + 



sinh (2 \q Q 


a) 






I 



(29) 

For $ 7^ 7r, the curves G(f2) and r g (£l) may become strongly asymmetric (as in Fig.4), 
the minimum in the gain curve is shifted away from Q = or multiple local minima may 
occurs; however the general rule of faster than v g \ group velocity near local minima of the 
power spectral gain curve holds. Figure 6 shows the behavior of the group delay at Q = 0, 
for the case of perfect phase reversal of the two grating sections (i.e., for $ = 7r), versus 
intensity of the pump wave. The figure clearly indicates the possibility of controlling the 
transit time of a spectrally-narrow signal pulse by the pump intensity. Notice that the group 
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delay becomes superluminal, i.e. smaller than L/cq, at I3 ~ 105 MW/cm 2 and negative at 
I 3 ~ 135 MW/cm 2 . To reach the latter regime, the pump intensity corresponds to an 
off-resonance power gain peak as large as 40 dB [see Fig.3(a)]. However, since the pulse 
spectrum is centered at Q = where G(Q) ~ 1, the main effect of the amplifier on pulse 
propagation near the gain dip is to advance the pulse in time, not to amplify it. A similar 
effect was predicted and observed for pulse propagation in an atomic gain-doublet amplifier 
|19| , |20|| . What physically happens is that in the fist section of the QPM grating the signal 
pulse is amplified, with the generation of the idler wave (003 — > u\ + cj 2 ); however owing to 
the phase reversal in the second QPM grating, a back conversion process (u>i + 002 — * W3) 
occurs in the second grating section. The result of such a cascading process is that the 
signal pulse is basically not amplified at the output of the crystal, however the phase delays 
suffered by its spectral components, as ruled by the causality condition [Eq.(24)], produce a 
temporal advancement with no appreciable pulse distortion. We have checked the prediction 
of pulse propagation based on the group delay analysis by direct numerical simulation of 
pulse propagation in a two-section QPM PPLN crystal with perfect phase reversal ($ = n) 
starting from Eqs.(8) in the limit of an undepleted pump. As a probing pulse, we assumed a 
transform-limited Gaussian pulse, tuned at Q — 0, with a pulse duration (FWHM) of 250 ps 
to ensure the spectrally-narrow pulse limit. Figure 7 shows the traces of the incident signal 
pulse at the entrance plane of the crystal z = (dotted line) and of transmitted pulses at 
the exit plane z = L (solid lines) for increasing values of the pump intensity ^3; the pulse 
intensity of the output waveforms are normalized to the peak intensity of incident pulse. It 
is remarkable that, at pump intensities corresponding to, e.g., curve 4, the transmitted pulse 
leaves the amplifier before the peak of Gaussian incident pulse has entered into the crystal. 
To make a quantitative estimate of pump power levels required to observe such effects, let us 
assume a Gaussian pump with a beam waist of ~ 200 //m; then curve 3 of Fig. 7 corresponds 
to a pump peak power of ~ 170 kW. Using a pulsed pump of duration (FWHM) of ~ 5 
ns, i.e. about twenty times longer than the probing pulses, a pump pulse energy of ~ 0.90 
mJ is required, which can be obtained using a frequency-doubled Q-switched Nd-based laser 
system as a pump source. 
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V. CONCLUSIONS 



In this paper we have presented a general analysis of the dispersive properties of QPM 
optical parametric amplifiers. In the limit of no pump depletion, the parametric interaction 
of idler and signal fields, that accounts for material dispersion at any order, may be described 
by coupled-mode equations which have a canonical form widely encountered in problems of 
inverse scattering and optical gratings and filters design. One of the main consequences 
of the analysis is that the spectral power gain curve of the amplifier defines univocally 
its dispersive curve through a Hilbert-like transform. This circumstance may be of major 
relevance in the design and synthesis of QPM amplifiers for applications in pulse shaping and 
control, and imposes unavoidable physical limits to the realization of dispersionless amplifiers 
with a flat spectral gain f3(|. Our analysis also demonstrates that engineered QPM grating 
profiles can simulate dispersive properties of resonant pulse propagation in inverted atomic 
media. In particular, in this work we have proposed and studied the possibility of observing 
and controlling superluminal and negative group velocities of picosecond optical pulses in 
a PPLN amplifier, which shows features similar to those found in atomic amplifiers with a 
gain doublet @ fE^ ^ . 



APPENDIX A: DERIVATION OF ENVELOPE EQUATIONS 

In this Appendix we derive the envelope equations [Eqs.(8a)-(8c)] given in the text by a 
multiple scale asymptotic analysis of the original equations [Eqs.(3a)-(3c)]. Such an analysis 
assumes that the nonlinearity of the medium is weak and that the spectral extent of the 
interacting fields is narrow enough such that in Eq.(5) the operator k 2 {ujQ + idt) — k 2 {uJo) can 
be treated as a perturbation term. With these assumptions, it is worth rewriting Eqs.(3a)- 
(3c) in the form: 



d 2 8 x 
dz 2 



d 2 8, 
dz 2 



^3 



+ k% = e | [k\ - k 2 {oo x + id t )] 81 - X i2) {^j ^ 2 *| (Al) 
+ k\8 2 = e I [k 2 2 - k 2 (u 2 + idt)] 82 - X (2) 2 ^ I (A2) 

+ k% = e I [kl - k 2 (u 3 + id t )] 8 3 - x {2) (^j 2 8 X 8 2 J (A3) 
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where e is a small parameter that defines the order of magnitude of the perturbation terms 
entering on the right hand side in Eqs.(Al)-(A3). The problem is to construct an asymptotic 
solution of the perturbed equations as e — > which is valid uniformly with respect to the 
spatial variable z. Therefore we look for a solution to Eqs.(Al)-(A3) in the form: 

S l = si 0) + eSi 1) + e 2 si 2) + ... (A4) 

(I = 1, 2, 3), and we require that the asymptotic expansion be uniformly valid. This condition 
can be satisfied by introducing multiple scales for z, i.e. by assuming that Si depends 
on Z , Z 1 , Z 2 , where Z = z, Z 1 = ez, Z 2 = e 2 z, ... The introduction of multiple 
spatial scales is fundamental to remove secular growing terms that arise in the perturbation 
expansion. Introducing expansion (A4) into Eqs.(Al)-(A3), using the derivative rule d\ = 
®z + e (^z ^i + 9 Zl d Zo ) + ••• an d collecting the terms of the same order in the equations 
so obtained, a hierarchy of equations for successive corrections to Si is obtained. At leading 
order, O(e ), one obtains: 

Q2p(0) 

-^r + *?*« (0) = °- (AS) 

(/ = 1, 2, 3). If we consider forward propagating waves, i.e. a copropagating interaction, the 
solutions to Eq.(A5) are given by: 

S^(Z , Z 1 ,...;t)=A l (Z 1 ,...;t) exp(ifc,Z ), (A6) 

where the amplitudes Ai depend on time and on slow spatial variables. At 0(e) one obtains: 

n2c(l) 

-^r + = g!" (AT) 

where the driving terms G^ in Eqs.(A7) are given by: 

G? = -2d Zo d Zl S^ + [kl - V + id t )} S^ - x < 2 > 2 Sfsf> (A8) 

G? = -2d Zo d z A° ] + H - k\u 2 + &)] 4 0) - X {2) (^) 2 W (A9) 
G? = -2d Zo d Zl St ] + [kl - k\u, + id t )] 4 0) - X (2) (£) 2 4 0) 4 0) . (MO) 
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To further proceed, let us assume that x^ ( z ) is a quasi-periodic function of z with a period 
A satisfying the QPM condition given by Eq.(7) in the text, and that A is much smaller 
than the crystal length and of the same order of magnitude as the wavelengths of interacting 
fields. It is then worth separating the fast and slow dependence of x^ 011 z by setting: 

oo 

X {2 \z)= xi 2) (^i)exp(-27rmVA) (AH) 

n=— oo 

(2) 

where the coefficients Xn of the Fourier series are allowed to vary on the slow spatial scale Z\ 
and satisfy the condition x~ n = Xn 1 * ■ The solutions to Eqs. (A7) are bounded with respect 
to Zq provided that the driving term does not contain terms oscillating like exp(ikiZ ). 
The solvability conditions allow one to derive the evolution equations of the envelopes A\ 
on the slow spatial scale Z\. Substitution of Eqs.(A6) and (All) into Eqs.(A8)-(AlO), 
after collecting the the terms oscillating like exp(ik[Z ), one obtains the following solvability 
conditions: 

-2^|| + [k\ - k 2 ^ + id t )] A, - x { 3 2 A*A 3 = (A12) 

~ 2 ^ 2 it + ^ 2 " e{uJ2 + idt) } M " X ^ (?) 2 AlM = (A13) 

~ 2lh Jz\ + ^ ~ e{uj3 + idt) ^ As " (?) 2 AlM = (A14) 

If we stop the asymptotic expansion at order ~ e, re-introducing the original spatial variable 
z instead of Z 1 in Eqs.(Al2)-(Al4) and setting e = 1, one finally obtains Eqs. (8) given in 
the text, where d e // is defined through Eq.(9). 



APPENDIX B: PROPAGATION OF SPECTRALLY-NARROW OPTICAL 
PULSES 

Let Ai(0, t) = exp(—iQ t)h (t) be the signal field envelope incident upon the amplifier 
with a frequency offset Qq from the carrier lo\. From Eq.(ll) the pulse waveform at the exit 
of the amplifier can be written as Ai(L,t) = hi{t) exp(—ifl t) , where: 

h L (t)=g(n + id t )h (t) (Bl) 
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and the operator on the right hand side in Eq.(Bl) is defined, as usual, by its power series 
expansion JO]]. We now assume that the spectral extent of the envelope ho(t) is narrow 
enough such that the gain curve g(fl) of the amplifier varies slowly over the pulse bandwidth. 
One can then expand g(fl) in power series around Q = Qq as follows: 



g(Q) = exp[\ng(fl)) = g exp 



1 ( d 2 \ng 
dQ 2 



(Q-n ) 2 + .. 



n 



g Bexp 



d\ng\ 



n 



where we have set: 



B(Q — Q ) = exp 



1 / d 2 \ng 



2 \ dtt 2 



(Q-n ) 2 + 



:b2) 



(B3) 



Substitution of Eq.(B2) into Eq.(Bl), after observing that the operator exp(— Ard t ) is equiv- 
alent to the time translation t — > t — Ar, one obtains: 



h L (t) = g B(id t )h (t - At) 



(B4) 



where we have introduced the "complex" group delay: 



Ar = -% 



dXwg 



dlnVG 
dQ 



n 



From Eq.(B3), the following asymptotic expansion holds for the operator Biidt 

1 / d 2 lna\ d 2 

B(id t ) = i-~ 1 M 



2 V dQ 2 J no dt 2 



(B5) 



(B6) 



where the dots involve higher-order time derivatives; substitution of Eq.(B6) into Eq.(B4) 
finally yields: 



h L {t) 



90ft0(t - AT) -|(^)^ (t - AT)+ ... (B7) 



At the leading order in the expansion, Eq.(B7) reduces to Eq.(25) given in the text. 
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Figure Captions 

Fig.l Schematic of parametric amplification of a weak signal wave at frequency uj\ 
in a QPM nonlinear crystal pumped by a strong CW pump at frequency u 3 (copropagating 
interaction) . 

Fig. 2 PPLN crystal geometry. In (a): the crystal comprises two uniform grating 
sections, each of length a, separated by a distance I with no grating structure. In each 
grating section, a +/- square periodic wave pattern of domains with period A is assumed. 
This yields d e ff = (2/71)^33 in the first grating and d e ff = exp(i<I>)(2/7r)<i33 in the second 
grating, where the phase $ is determined by the relative phase shift of square waves in 
the two grating sections. In (b) it is shown the QPM square wave profile x^K z ) that 
corresponds to $ = it. 

Fig. 3 Behavior of spectral power gain (a) and group delay (b) for $ = it. The 
pump intensity is ^3 = 135 MW/cm 2 . The other parameter values are given in the text. 
Solid and dashed curves refer to the results obtained by use of Sellmeier equations, i.e. 
taking into account dispersion at any order, and by use of approximate equations (22) given 
in the text, respectively. The inset in (b) shows an enlargement of the group delay near 
resonance. 

Fig.4 Same as Fig.3, but for $ = it/2. 
Fig. 5 Same as Fig.3, but for $ = 0. 

Fig. 6 Behavior of group delay r g versus pump intensity I3 at the center of the am- 
plifier gain dip. Parameter values are the same as in Fig.3. 

Fig. 7 Traces of the intensity of transmitted signal pulse (solid curves) at the exit 
plane z = L of the amplifier for a few values of pump intensity. Curve 1: ^3 = 0; curve 2: 
J 3 = 108 MW/cm 2 ; curve 3: I 3 = 135 MW/cm 2 ; curve 4: J 3 = 162 MW/cm 2 . The dashed 
curve is the trace of incident Gaussian pulse at the input plane z = 0. 
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